The precision of slow-roll predictions for the CMBR anisotropics 
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Inflationary predictions for the anisotropy of the cosmic microwave background radiation (CMBR) 
are often based on the slow-roll approximation. We study the precision with which the multipole 
moments of the temperature two-point correlation function can be predicted by means of the slow- 
roll approximation. We ask whether this precision is good enough for the forthcoming high precision 
observations by means of the MAP and Planck satellites. The error in the multipole moments due to 
the slow-roll approximation is demonstrated to be bigger than the error in the power spectrum. For 
power-law inflation with ns = 0.9 the error from the leading order slow-roll approximation is « 5% 
for the amplitudes and ~ 20% for the quadrupoles. For the next-to-leading order the errors are 
within a few percent. The errors increase with |ns — 1|. To obtain a precision of 1% it is necessary, 
but in general not sufficient, to use the next-to-leading order. In the case of power-law inflation this 
precision is obtained for the spectral indices if |ns — 1| < 0.02 and for the quadrupoles if jris — 1| < 0.15 
only. The errors in the higher multipoles are even larger than those for the quadrupole, e.g. ~ 15% 
for I — 100, with ns = 0.9 at the next-to-leading order. We find that the accuracy of the slow-roll 
approximation may be improved by shifting the pivot scale of the primordial spectrum (the scale 
at which the slow-roll parameters are fixed) into the regime of acoustic oscillations. Nevertheless, 
the slow-roll approximation cannot be improved beyond the next-to-leading order in the slow-roll 
parameters. 
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I. INTRODUCTION 

High quality measurements of the Cosmic Microwave 
Background Radiation (CMBR) anisotropics have been 
pubhshed recently by the balloon borne experiments 
BOOMERanC § and MAXIMA-1 g. A large num- 
ber of multipoles (26 < / < 625 for BOOMERanC and 
36 < Z < 785 for MAXIMA-1) has been covered by both 
experiments. During the next years, high precision mea- 
surements will be performed by the MAP and Planck 
satellites Inflation jj] provides a mechanism to pro- 
duce the primordial fluctuations of space-time and mat- 
ter H"-^, which lead to the CMBR anisotropics and to 
the large scale structure. This mechanism rests on the 
principles of General Relativity and Quantum Field The- 
ory. It thus can be expected to get a hand on the physics 
of the very early Universe with help of the upcoming high 
precision measurements. 

The CMBR anisotropics are most conveniently ex- 
pressed by the multipole moments C/. The compu- 
tation of the multipole moments requires the knowl- 
edge of the primordial spectrum and the transfer func- 
tions. The latter depend on the cosmological parameters 
Hq, r^M, ^^A, • ■ •■ The transfer function characterizes the 



evolution of cosmological perturbations during the radi- 
ation and matter epochs. The primordial spectrum is 
predicted by inflation and depends on the evolution of 
the long wavelength perturbations during inflation and 
reheating. It can be predicted from a given model of 
inflation. 

In this article, we will restrict our considerations to 
slow-roll inflation with one scalar field. This represents 
only a first step towards a more general study. Our aim 
is to address the following problems: What is the preci- 
sion of the predicted multipole moments from the slow- 
roll approximation? Is this precision sufficient to reach 
the level of accuracy expected from the planned observa- 
tions? Can the slow-roll approximation be improved to 
arbitrary precision? 

So far, the precision of the predicted power spectrum 
has been examined by Grivell and Liddle However, 
the power spectrum is not directly observable whereas 
the C;'s are. We show that the error from the slow-roll 
approximation is important in the multipole moments. It 
is bigger than the error in the power spectrum. It turns 
out that the next-to-leading order slow-roll approxima- 
tion ||l^ is compulsory, but it may not be sufficient to 
reach an accuracy of a few % or less. 
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Wang, Mukhanov and Steinhardt |T|] have shown that 
predictions based on the time delay argument Q or on 
the horizon crossing/Bessel function approach (i.e. the 
slow-roll approximation) ||l^,|l2| are not reliable for gen- 
eral models of inflation. There have been various at- 
tempts to improve the slow-roll approximation to higher 
orders, see e.g. JTo| , p^ , p^p^ . The conclusion of Wang et 
al. [|ri[ has been contested by Copeland et al. "We 
. . . conclude that any theoretical errors from the use of 
the slow-roll equations are likely to be subdominant" . We 
show in this work that this claim is not correct unless the 
slow-roll parameters are extremely small. Typically, we 
find that the slow-roll parameters must be less than 0.01 
in order for the next-to-leading order to reach the level 
of precision of MAP or Planck. This means that there 
are models where the slow-roll error is dominant and the 
slow-roll approximation is valid. We find in agreement 
with the analysis in that a slow-roll approximation 
that goes beyond the next-to-leading order cannot ex- 
ist. All higher order corrections are thus meaningless. In 
the derivation of this result we close a gap in the proof 
of the next-to-leading order equations. For some reason 
this gap was not noticed before in the literature. For 
this purpose we use and generalize a new family of exact 
solutions, which was recently found by Starobinsky . 

The scope of this paper is to quantify the error from 
the slow-roll approximation. We compute the scalar mul- 
tipole moments and the ratio R = Cj/Cf for power-law 
inflation for which the exact result is known. Then, we 
calculate the same quantities for the same model but in 
the context of the slow-roll approximation. The compar- 
ison of the two results provides an estimate of the error 
made by using the slow-roll approximation. We do not 
convolute this error with the uncertainties in the trans- 
fer functions. For the sake of clarity and simplicity we 
only make use of the transfer functions in the long wave- 
length limit. This approximation only mildly affects the 
estimates of the error in the multipolc moments. Then, 
we compare the slow-roll errors at leading and next-to- 
leading order to the cosmic variance. Binning several 
multipoles together allows to reduce the cosmic variance, 
but does not reduce the slow-roll error. We find that the 
slow-roll error is hidden in the cosmic variance only for 
very small values of the slow-roll parameters (< 10"^). 
We propose to reduce the slow-roll error by optimizing 
the pivot scale (the scale at which the slow-roll param- 
eters are fixed) of the spectrum. However, this method 
is not sufficient to hide entirely the slow-roll error in the 
cosmic variance. 

This article is organized as follows: in the next section, 
the theory of cosmological perturbations and the calcula- 
tions of the CMBR anisotropics are reviewed. Then, the 
low-^ multipole moments are computed exactly for power- 
law infiation (Sec. Ill) and approximately for slow-roll 
infiation (Sec. IV). Comparison of the two results allows 
us to test the precision of the CMBR multipoles obtained 
from the slow-roll approximation in the last section. The 
slow-roll errors are shown to be observationally signifi- 



cant by comparing them with the cosmic variance. We 
set c = h = 1 throughout the paper. 



II. FROM QUANTUM FLUCTUATIONS TO 
CMBR ANISOTROPIES 

The line element for the spatially flat Friedmann- 
Lemaitre-Robertson- Walker background plus perturba- 
tions can be written as 1^: 

ds^ = a^{rj){-{l - 2(l>)dT]^ + 2{d^B)dx'dri + [(1 - 2i;)6,j 
+2didjE + hij]dx'dx^} . (1) 

In this equation, the functions 4>^ B, ip and E repre- 
sent the scalar sector whereas the tensor /ly , satisfy- 
ing hi^ — hij'-' — 0, represents the gravitational waves. 
There are no vector perturbations because a single scalar 
field cannot seed rotational perturbations. The confor- 
mal time rj is related to the cosmic time thy dt — a{iiii)drj. 
It is convenient to introduce the background quantity 
7(77) defined by 7 = —H/H^, where a dot means differ- 
entiation with respect to cosmic time and H is the Hub- 
ble rate, H = a/a. Using conformal time we may write 
7 = 1 — H' /Ti,^, where Ti = a' /a, and a prime denotes 
differentiation with respect to the conformal time. 

We assume that inflation is driven by a single scalar 
field. For the perturbations we introduce gauge-invariant 
variables PTUq], which reduce the equations of motion, in 
the small scale limit, to equations of harmonic oscillators 
|]l8|-^^. In the tensor sector (which is gauge invariant) 
we define the quantity /it for each mode k according 
to hij = {^T / a)Qij{k), where Qij{k) are the (transverse 
and traceless) eigentensors of the Laplace operator on the 
spacelike sections and fc^ is the corresponding eigenvalue. 
Gravitational waves do not couple to scalar fields. Thus 
the equation of motion is given by [Esl: 



a 
a 



Mt = 



(2) 



The scalar sector is gauge dependent and the scalar per- 
turbations of the metric are coupled to the perturbations 
of the stress tensor describing the matter. Fluctuations 
in the stress tensor involve perturbations in the energy 
density, Sp, and in the four velocity, Su^ = {—(/)/ a, v^/ a). 
We describe perturbations in the density contrast by the 
gauge invariant quantity 6 = Sp/p+ {p' / p){B — E'). Per- 
turbations in the velocity can be written as Vi = dtw+Wi. 
Since we are interested in the scalar sector, only the 
first term has to be taken into account. We choose to 
work with the gauge invariant quantity v = w + E'. 
Scalar perturbations of the geometry can be character- 
ized by the two gauge invariant Bardeen potentials = 
(t>+{l/a)[iB-E')a]' and *g = -U^B - E') 0, where 
Q{k) is a scalar harmonic. During inflation, the Universe 
is dominated by the scalar field (/? = (/3o(f?) + 'Piiv)Q- 
Fluctuations in the scalar field are characterized by the 
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gauge invariant quantity 5f = 'pi+ ip'q{B — E'). In this 
simple case, the time evolution of fluctuations can be 
reduced to the study of the equation of motion for the 
variable /is = —^/2Ka[S(f + {(p'q/H)^], where k = 8ttG. 
Its equation of motion is very similar to that of the grav- 
itational waves l[13,E0i: 



(«V7) 



MS = 



(3) 



The integration of (||) and ^ leads to the primordial 
spectrum of the fluctuations. For the initial conditions 
we assume that the scalar and tensor perturbations are 
in the quantum vacuum state when the scale of interest 
was well within the Hubble radius (l/Zcph <^ c/H) during 
the early stages of inflation. Therefore all fluctuation 
variables are quantum operators during inflation. After 
inflation, the Universe is filled with baryons, photons, 
neutrinos and (cold) dark matter. For that epoch, the 
perturbed Einstein equations cannot be reduced to the 
simple form of Eqs. (||) and and need to be integrated 
numerically. This leads to the transfer functions. 

The cosmological perturbations induce anisotropics in 
the temperature of the CMBR, which have been detected 
by COBE H first. This is the Sachs- Wolfe efi'ect [||. 
Since it does not depend on the photon frequency, the 
black body shape of the photon spectrum is conserved 
from the last scattering surface to its observation today 
p^ . The measured anisotropics in the photon intensity 
translate into anisotropics in the temperature of the black 
body. 

For the temperature fiuctuations we introduce the ab- 
breviation A(e) = {ST/T){e), where characterizes the 
direction of the beam on the celestial sphere. The contri- 
butions of the scalar and tensor perturbations are given 
by: 



1 f" 8 

A^{e} = -5^ + $ - e'd^v + / dfj—{^ + *) , 



(4) 



(5) 



The first three terms of the scalar contribution are evalu- 
ated on the last scattering surface, i.e. at ?7iss. They repre- 
sent the intrinsic fluctuations, the Sachs- Wolfe effect and 
the Doppler effect. The forth term is the so-called inte- 
grated Sachs- Wolfe effect. The integration is performed 
along the photon trajectory, which is parameterized by 
the conformal time here, rjo denotes the conformal time 
at observation today. 6^ is the perturbed density con- 
trast of the photons and v the perturbed velocity of the 
photon fluid. For large angular scales only the flrst two 
terms are important. For isentropic (sometimes called 
adiabatic) perturbations the scalar part reduces to: 



Im + {■■■) 



(6) 



Usually, the CMBR anisotropics are expressed through 
the multipole moments C; . The C; are the coefficients in 



an expansion over Legendre polynomials of the CMBR 
temperature two-point correlation: 

(A«-T(ei)AS^T(.)^ = i-^(2Z + l)Cf'^P,(cos5) , (7) 



where cosi5 = ei • 6*2. The brackets () denote the averag- 
ing over many ensembles. Averages over many ensembles 
cannot be replaced by spatial averages on the celestial 
sphere due to the lack of ergodicity of the stochastic pro- 
cess A(e), see Ref. 25|. If, nevertheless, we do this the 
error made can be quantified by means of the cosmic 
variance. 

The computation of the multipoles for a given model 
requires the knowledge of the initial spectrum of the fiuc- 
tuations and of the transfer function. The power spec- 
trum of the Bardeen potential is defined in terms of the 
two-point correlator for the operator $(77, x): 



f°° c\k sin kr 

(0|$(7y,x)<f(r,,x + r)|0) ^ ^t_k^P^^rj,k) . (8) 

Similarly, the power spectrum of gravitational waves is 
defined as: 



|/i,,(77,x)/^*^"(77,x + r)|0) 



dfc sin kr 
k kr 



(9) 



A priori, the primordial power spectra are time depen- 
dent quantities. However, for the multipoles between 
I = 2 and I = 2000, we are interested in scales which 
are well beyond the horizon at the end of inflation. In 
a first approximation for those scales the power spectra 
do not evolve in time during infiation and they can be 
written as: 



fc3p<j,(fc) = A^(fco) ( ^ 
ePh{k)^M^{ko) 



k 

^0, 



(10) 
(11) 



where the spectral indices ns, and the amplitudes 
v4g, A\. are independent quantities and fco is an arbi- 
trarily fixed scale which is introduced to link various 
notations in the literature. The spectral indices can 
also be determined from rig — 1 = dln(fc'^P$)/dlnfc and 
TT-T = dln(fc3p^)/dlnfc. 

An accurate calculation of the multipole moments re- 
quires numerical computations. However, for small /, the 
approximate equation (^ can be used. For density per- 
turbations Eql this leads to: 



An 



dfc 



jf{kriss)T^{kriss ~^ 0) 



xAl ( — 

fco 



ns-l 



(12) 
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where ji is the spherical Bessel function of order I and 

rjP^ = a{r]o)riss = ao(??o - ??iss) « ao?7o ~ 2i?H is the 
comoving hne-of-sight distance to the last scattering sur- 
face. 

T$(fcriss — > 0) is approximately the transfer function 
for superhorizon modes. It is k independent and there- 
fore only the amplitude is modified but not the spectral 
index. The domain of validity of the latter approximation 
can be evaluated as follows. In the integral (12) the main 
contribution comes from the modes around kriss ^ I + 1. 
We use that to estimate for which multipole moments 
the k independent transfer function is good enough. The 
mode whose wavelength is equal to the Hubble radius 
today, i.e. such that 27rao/fc = Ih{,Vo)j has krjo — 47r. 
Therefore the constant superhorizon transfer function is 
a reasonably good approximation if [l + l)/nss ^ 47r/77o, 
that is to say I <C 10. This is a very optimistic estimate 
since it does not take into account the first approxima- 
tion, made in Eq. (I), on which the vahdity of Eq. (|l|) 
rests. 

With the above approximations the low-/ multipoles 
can be calculated exactly pffl. The result reads: 



_3/, r[(3-ns)/2]r[/ + (ns-l)/2] 



r[(4 

x(fconss)'""^ 



ns)/2]r[/- 
9 ' 



("s - l)/2] 



(13) 



where Ag = AgT,^{kriss — > 0). For gravitational waves 
we obtain the following expression p7[ : 



9tt 



il-i)lil + i)il + 2) (konssY 



^|/K2/)|'At2/"-, 

y 



where the function /; (y) , y 



is defined by: 



Jo x(y - xY 



(14) 



(15) 



The superhorizon transfer function for gravitational 
waves does not appear explicitly because it is equal to 
one. As a consequence we can write Alj, = At^. The com- 
putation of C]' is more complicated than the calculation 
of Cf . The integral /; can be calculated exactly in terms 
of special functions, see Ref. However, the second 

integration over k cannot be performed analytically and 
we must rely on numerical integration. 

Below we will be interested in the ratio of tensor to 
scalar quadrupole contributions [p8|-pl|: 



(16) 



Expressed in terms of the tensor spectral index this is 
the so-called consistency equation of inflation. 

We have seen that the calculation of Cf requires the 
knowledge of the transfer function and of the primordial 



spectrum. In principle, T,^(kr\ss) is known accuratel y a s 
the result of numerical calculations, e.g. see Refs. [p2[ . 
When we calculate the multipoles using Eqs. ( p^ ) or 
( |l^ ) we make two approximations: a long wavelength 
approximation for the transfer function and we neglect 
the contribution of radiation (pure matter assumption) 
to the expansion of the Universe at photon decoupling. 
The long wavelength approximation results in neglecting 
other contributions besides (^) in the Sachs- Wolfe efllect 
for scalars and in considering that the tensor and scalar 
superhorizon transfer functions are constant. The pure 
matter assumption results in a small error in the numer- 
ical value of r<i>(/criss ^0). However, for small values of 
I these errors are small. 
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FIG. 1. Error due to the long wavelength approximation 
in the transfer function for the scalar multipoles with a flat 
primordial spectrum. The exact multipoles are calculated 
by means of the CMBFAST code and are normalized to the 
quadrupole. 

In order to test this claim quantitatively and to quan- 
tify the contribution to the total error coming from the 
transfer function (see Fig. |^), we compute the scalar 
multipole moments for low / numerically with CMB- 
FAST for the following values of the cosmologi- 
cal parameters: ~ 50km/s/Mpc, 57o = l,ilcDM — 
0.95, r^B = 0.05. We compare them to the multipole 
moments given by Eq. ( p^ with a constant transfer 
function. The code CMBFAST automatically normal- 
izes to the COBE result [ p2| . The result is expressed 
by the band powers {STJT^ = l{l + l)C;/(27r), where 
To w 2.73 K is the average temperature of the CMBR. 
For a flat (ng = 1) primordial spectrum, CMBFAST 
gives 5T2 ~ 27.5fj,K or Qrms-ps ~ 17.8fiK, where 
the quadrupole rnis fluctuation is given by Qrms-ps = 
Tq\/ (5/47r)C2. We normalize the amplitude As in 
Eq. (n2|) is to the latter value of the quadrupole. In 
Fig. ^pwe plot the differences of both calculations, di- 
vided by the CMBFAST results, and express this num- 
ber as the error in %. The error in the quadrupole 
(O) vanishes "by construction". Equation (03) shows 
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that STi = ToVA^/3(fconss)(^-"^^/2 = Orms-PS^WS, 

whereas the CMBFAST-JT/ is ^-dependent, despite both 
band powers are calculated from the same primordial 
spectrum. The difference between both band powers 
is exclusively due to the use of different transfer func- 
tions and to the neglection of the Doppler and integrated 
Sachs- Wolfe effects. In this way, we can isolate and esti- 
mate the error coming from the long wavelength approx- 
imation, being given that the spectrum is normalized to 
COBE. 

A similar study has been done in Rcf . |Q . The errors 
given in that article differ from those obtained here be- 
cause a different normalization is used. In Ref. |3^, the 
spectrum is normalized to the multipole moment Cio in- 
stead of the quadrupole. As a consequence, in that case 
the error in Cio vanishes "by construction" . 

Figure |^ confirms the importance of the transfer func- 
tion and the analytical estimates made at the beginning 
of this article. The error is below 1% only for I < 4. For 
Cio, which is often used to normalize the spectrum, the 
effect of the subleading terms in k is already 5%. The 
error from the pure matter assumption has not been fully 
accounted for by this method, because we do not test the 
error in the numerical value of r$(A:riss —> 0) when we 
normalize the quadrupole to the COBE result. Since this 
error is a pure overall numerical factor, it does not affect 
our conclusions. 



III. PREDICTIONS OF POWER-LAW 
INFLATION 

In this section, we turn to the study of power-law in- 
flation. This model is of particular importance because 
it allows to calculate all quantities of interest exactly. 
Moreover, this exact result is at the basis of the slow-roll 
approximation. 

Power-law inflation is given by the following solution 
for the scale factor and the scalar field: 



^^^, + !I^Jlil+p)lnM, (17) 



where mpi is the Planck mass and is the initial value 
of the scalar field at conformal time r]i. In this model 
inflation occurs if /3 < — 2 (we do not consider the case 
where — 2 < /3 < — 1 which cannot be realized with a 
single scalar field) . The quantity Iq has the dimension of a 
length and its value will roughly determine the amplitude 
of the CMBR fluctuations today. In the particular case 
of power-law inflation, the function jij]) is a constant 
equal to {2 + + (3). For -oo < P < -2, 7 goes from 
one to zero, this last value corresponding to the de Sitter 
spacetime. The scale factor and scalar field of Eqs. ( |l7| ) 
are solutions of the Einstein equations for the scalar field 
potential: 



V{(p) = Vi exp 



mpi 



Vii'P - Vi) 



(18) 



where Vi is the value of the potential at 77;. 



A. Density perturbations 

The effective potential for density perturbations. Us 
(«V7)'7(aV7), see Eq. (|), reads: 



Usiv) 



(19) 



This simple form of the potential allows an exact inte- 
gration of Eq. (^). The solution is expressed in terms 
of Bessel functions. This provides the initial power spec- 
trum, i.e. Ag and ng. In order to evolve the superhorizon 



spectrum, we can rely on the conservation law 21 
for the quantity: C = -f <i>)/7 -|- <1>. This gives 

the superhorizon transfer function: T$(fcriss ^ 0) = 
[9(2/3 -f 3)2]/[2572(l + I3f]. Then the amplitude of the 
scalar quadrupole and the spectral index take the form: 



As(fco) 



where 



9 



257r7 



Tls-1 



ng = 2/3 + 5 



r(-/3-i/2) 

2/5+1 



1 - 37 

(20) 



(21) 



which is unity for /3 = — 2. As expected, the amplitude of 
scalar perturbations is roughly determined by the ratio 
Ipi/lo- Very often the final spectrum is expressed in terms 
of the Hubble rate at some time 77*, instead of the scale ^q. 
We have = H{rj^) = -[(l-h/3)//o]h*r^"^. Therefore 
the amplitude reads: 



9 



257r7 



/(/3)|l+/3|2(''+i) 



fcn 



(22) 

The amplitude Ag is displayed as a function of 7 in Fig. |^. 
It diverges in the de Sitter limit 7 — + 0. COBE measured 
the spectral index to be ns = 1.2 ± 0.3 The l(2)cr 
value ns — 0.9(0.6) corresponds to 7 « 0.048(0.167). 



B. Gravitational waves 

The calculation of the spectrum for gravitational waves 
is performed along the same lines as above. The effective 
potential is the same as for density perturbations, i.e. 
Ut = a" /a = /3(1 + (3)/rf' ■ Since the superhorizon trans- 
fer function is equal to one, A'y and rtx can be written 
as: 



A-rikQ 



72 16 

In TT 



riT = 2/3 + 4 = - 



27 

1-7 



(23) 
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For power-law inflation, the relation ng = ut + 1 holds. 
In terms of the amplitude is given by 



r2(3:^)r(4- 



(24) 



Figure || shows the scalar and tensor amplitudes (^2|) and 
(H), respectively. For 7 > 9/400 = 0.0225 the tensor 
mode dominates. 
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FIG. 2. The amplitudes of scalar and tensor perturba- 
tions. In the de Sitter limit 7 ~+ the scalar amplitude di- 
verges. For larger values of 7 the perturbations are dominated 
by the tensor mode. 



C. Multipole moments 

The multipole moments predicted by power-law infla- 
tion can easily be computed from Eqs. ( p^ and The 
quadrupoles are displayed in Fig. ||. Compared to the 
amplitudes the importance of the tensor mode is slightly 
suppressed, it becomes the dominant mode at 7 > 0.07, 
which corresponds to ns ^ 0.85. 



Power-law inflation 




0.05 0.1 0.15 0.2 

7 

FIG. 3. The quadrupole moments of scalar and tensor 
perturbations. 

We calculate the ratio R for power-law inflation: 
R = 13.867F[nT(7)] = -6.93-^^i^(nT) , (25) 
where the function _F(nT) is given by: 



F{nT) = 496.1 X 2 



1 — /It 



2 ) 



r(2 - nT)F(2 + 



dkk 



riT-l 



WW 



(26) 



In this expression we have used the equation ng = nT + 1, 
valid for power-law inflation only, to express everything 
in terms of nx. We have dkk-^lhik)]^ = 2.139x10"'^ 
such that F{nT = 0) = 1. Notice that the factors fcoHss 
and kQ/{aifH^) cancel in R because ng = nx + 1. R 
versus 7 is plotted in Fig. |4[ This plot demonstrates 
that within the 2a error bars of COBE, there is a large 
parameter space where the tensor mode dominates the 
scalar modes, see e.g. Refs. |2^,0 for a more detailed 
discussion. 




FIG. 4. The tensor to scalar ratio of the quadrupole mo- 
ments. 



IV. PREDICTIONS OF SLOW-ROLL INFLATION 

For a general model of inflation exact solutions are not 
available. Generically, the potentials Us and Ut are dif- 
ferent but nevertheless their shape is similar. A sketch of 
the generic form of J7g and is displayed in Fig. ||. The 
details of the realistic reheating transition are not taken 
into account in this simple figure. During the radiation 
dominated era the potential goes to zero, since a oc rj. 

For a given mode fc, the inflationary epoch can be di- 
vided into three stages, see Fig. |[ In region I the mode k 
is subhorizon. In that case the effective potential is small 
compared to fc^. In the limit k/{aH) — s- 00 for fixed k, 
the vacuum fluctuations are given by, see Ref. |21|: 



-ik{ri — rii) 



Pl- 



2k 



(27) 



respectively. In region III the mode is superhorizon. In 
the hmit k/{aH) ^ at fixed fc, the potential term is 
dominant, and the "exact" solutions read: 



MS ('7) = Cs{a^){'q) X 

/"? 1 



{a^l){'n) 



{a'^j){fj) df] dry 



(28) 
(29) 
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Usually, density perturbations are described in terms of 
the Bardeen potential $ instead in terms of fis- The 
order fc^ term is necessary to obtain the leading or- 
der expression for the Bardeen potential, since <i> = 
[H-//{2P)][fj.s/{ay/^)]', see Refs. Thus, in region 

III, the superhorizon Bardeen potential is given by: 



$(?7) 



CsH 
2a2 



a^7 dr] 



(30) 



Our aim is to calculate the spectra at the end of inflation, 
i.e. in region III. The time dependence of the solutions 
in this region is known and the difficulty lies in the cal- 
culation of the constants Cs and Ct . Since the solutions 
are uniquely determined in region I, this amounts to join 
the super- and subhorizon solutions. Therefore we need 
to know the behavior of the perturbations in region II. 




1 2 



S) 



(34) 



The slow-roll conditions are satisfied if e and 5 are much 
smaller than one and if ^ = 0{e^ ,5"^ ^e5). From Eqs. 
(p^), it is clear that this amounts to consider e and 5 as 
constants. This property is crucial for the calculation of 
the perturbations. 

For power-law inflation the slow-roll parameters sat- 
isfy: 



= (5< 1 , e = 



(35) 



Therefore the slow- roll conditions are fulfilled if e ^ 1, 
that is to say if (3 is close to —2 (scale invariancc). In 
fact, the slow-roll approximation is an expansion around 
power-law inflation with < —{(3 + 2) <^ 1. To illustrate 
this point, let us consider the exact equation: 



1 



da- 



(36) 



If we assume that e is a constant, the previous equation 
reduces to aH « — (l-|-e)/?7. This is equivalent to a scale 
factor which behaves like: 



a(ry) w lo\-q\ 



-l-e 



(37) 



Interestingly enough, the effective power index at leading 
order depends on e only. 



A. Density perturbations 



FIG. 5. Sketch of the effective potential for density per- 
turbations and/or gravitational waves during inflation and 
radiation. 

A popular approach is the slow-roll approximation 
. The idea is that there was an epoch during infla- 
tion where the scalar field was rolling down its potential 
V{^p) very slowly. Under certain conditions (see below) 
this is close to the behavior during power-law infiation 
and the exact solutions from power-law inflation are used 
in region II to interpolate between the sub- and super- 
horizon solutions. 

Slow roll is controlled by the three (leading) slow-roll 
parameters (see e.g. Ref. Q) defined by: 



5 = 



3^f^ 

Jf_ ^ _ 
Hip 

-6 



V 



H 



2He 



H 



(31) 
(32) 

(33) 



We see in particular that 7(77) = e in region II. The equa- 
tions of motion for e and 5 can be written as: 



The effective potential of density perturbations can be 
calculated in terms of the slow-roll parameters exactly. 
The result is: 



Usiv) - «'^'[2 - e + (e - (5)(3 - J) + e] 



(38) 



-2(1 + 2e) and 
6e-3(5)?7-2. 



In the slow-roll approximation o?H^ rj 
the effective potential reduces to [/s ~ (2 
Since e and 5 must be seen as constants in the slow-roll 
approximation, the equation of motion (^) is of the same 
type as in power-law inflation. The solutions are given 
by Bessel functions: 



/is - (fcr/)i/2[i?iJ (3o(fcr/) -f B2J_„(=.)(M], 



whose order is given by 



,(sr) 



2€ + 5 . 



(39) 



(40) 



A comment is in order here: The potential depends 
on the scale factor and its derivatives only. One could 
think, looking at Eq. (p7|), that Us also depends on e 
only. This is not the case. The reason is that Us contains 
terms like e/e (for example) which are linear in 5, see 
Eqs. (^3). First one must calculate all derivatives, replace 
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them with their expression in terms of e and S, and only 
then consider that the slow-roll parameters are constant. 

We would also like to stress that keeping higher orders 
in e does not make sense. If terms of quadratic order in 
the slow-roll parameters are kept, the solution for density 
perturbations in region II can no longer be expressed in 
terms of Bessel functions. This is because the slow-roll 
parameters can no longer be considered as constant in 
time, see Eqs. (|3^). Therefore any considerations at this 
order in the framework of the slow-roll approximation is 
meaningless. The same conclusion has been obtained by 



Wang, Mukhanov, and Steinhardt |11 



Let us now calculate the constant Cg ■ The first step is 
to match the solutions of region I and II. This procedure 
fixes Bi and B2. Using Eqs. ( p7|) and (|2^), one obtains: 



and Bi 

>r)\ 



sm Trv, 



= 2z7rZpi exp[wg^'''' (7r/2) — 
Note that Bi and B2 do 



B1/B2 = -e^^'^s 

i(7r/4) -I- ikr]i]/{Vk 
not depend on the time at which the matching between 
regions I and II is performed. The joining between re- 
gions II and III remains to be performed at some time 
T^g, which will be fixed below. Expanding everything up 
to next-to-leading order in the slow-roll parameters, one 
obtains: 



\Cs\ 



1 - 2(C + ln/c)(2e - 6) 
+ 2(J-e)ln|77s|lr 



(41) 



with CEE7£;-hln2-2« -0.7296, 7e « 0.5772 be- 
ing the_Euler constant. The Bardeen potential given in 
is now completely specified. Note that 7(7/) in 
is a time dependent function, evaluated in re- 
gion III, whereas e in Eq. (41) is a constant parameter, 
which is fixed by e = 7('7s)- 

For scalar perturbations it is useful to evaluate the 
quantity 




c 



MS 



(42) 



which is a constant for the dominant mode at superhori- 
zon scales |^,||,^. The quantity — C is denoted TZ in 
Ref. [|l2|. Instead of expressing the spectrum in terms 
of the ratio l-pi/lo, it is usual to write it in terms of the 
Hubble rate at some time 77*. Of course, there is nothing 
deep in this choice and one could have kept working with 
Zpi/Zo- A priori, the value of 77* is arbitrary and could ei- 
ther be in regions I, II or III. However, in order to make 
contact with the literature, we will assume that is in 
region II. Then, in the slow-roll approximation, the value 
of H[rj^) can be written as: 



EE H{r]., 



[1 + e(l + In \r].. 



(43) 



In Ref. |lj], ?7* is the time which satisfies the relation 
a(77,)i?(?7,) = k for each mode k. In other words, we have 



?7* = ri^(k). In this article, we adopt another convention 
and choose 77^ such that it is not a function of k. Then, 
a straightforward calculation gives 



k^P^{k) = 



1 - 2e 



2{5-e) In 



2[C + lnfc|77,|] (2e-(5) + 

(44) 



The matching time 773 remains to be fixed by a physical 
argument. To our knowledge, this issue has been over- 
looked in the literature so far. All works on the slow-roll 
approximation, starting with Ref. [|lO|, have tacitly as- 
sumed that 773/77* = 1, without further justification. A 
priori, an equally good choice would be, for example, 
when the mode a*ff* crosses the effective potential, i.e. 
when (1 -f 2e)/77^ = Us{r]s)- It is easy to show that this 
boils down to the choice rj^/rj^, = \J2. It is important to 
realize that different choices for the ratio rj^/rj^: lead to 
different observational predictions. Although a change in 
?7s would not change the spectral index, it would change 
the amplitude of scalar perturbations and the ratio of 
tensor to scalar contributions R. 

The missing physical argument comes from a new fam- 
ily of exact solutions which has a slow-roll regime in a 
certain limit. One exact solution is of course power-law 
inflation, but it does not help for the purpose of fixing 
773/77*, because the spectrum does not depend on 773/77* 
for 5 = e. These solutions are found by the ansatz 



r ^ 

m 



(45) 



where A and a are two free parameters. This defines 
a two-parameters family of exact solution. Note that 
this family is not equivalent to power-law inflation. The 
power-law model {A = ^qvtI/^), a = —1 — /3) is a sub- 
class of this two-parameter family. Of course, this is be- 
cause 0(77) (X |77|^" is just a solution of Eq. (p5|), viewed 
as a second order differential equation for the scale factor, 
but not the general solution. The limit A to zero and a 
close to one gives a slow-roll inflation model. The particu- 
lar case a = 1 was already found recently by Starobinsky 
]l6[ . This one parameter family of solutions is character- 
ized by a flat spectrum, 715 = 1. Eq. (45) is a generaliza- 
tion of Starobinsky 's original ansatz |16|. The spectrum 
may be calculated exactly to read 



k^P^{k) 



(46) 



The case a = 1 gives k^PAk) — /|j/(7ryl^) and coincides 
with the result of Ref. [|||. 

We now need to calculate the slow-roll spectrum for 
this new class of solutions. A comparison with Eq. ( |4^ ) 
will allow us to fix the ratio t/s /77* . Let us first determine 
the slow-roll parameters. In the slow-roll approximation 
we find 
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V 



-(l + 2e-5) , 



(47) 



whereas insertion of the ansatz (|45|) into this equa- 
tion gives (0-^7) '/(a^/7) = —a/rj. Therefore, one has 
2e = {a — 1) + 5 and especiaUy 2e = (5 if a = 1. It is 
interesting to note that we no longer have the relation 
e — S typical of power-law inflation. Let us also em- 
phasize that the two-parameter family is the only fam- 
ily of exact solutions which permits a slow-roll approx- 
imation. Eq. (47) is a necessary condition for the va- 
lidity of the slow-roll approximation. This equation can 
be viewed as a first-order differential equation for the 
quantity a^/j. Integration of this equation leads to the 
ansatz given in Eq. (|4^). Therefore our determination 
of the ratio rys/??* is general. The value of A in the 
slow-roll limit is obtained from A = Q^/7|7?|" and is ex- 
pressed in terms of H^, with the help of (pq). This gives 
A"^ = ei/-2[l-h2e-h2(2e-(5)lnj77,|]. Thus we obtain the 
slow roll spectrum from Eq. (Eq): 



7re 
- 2(C 



[1 - 2e 

\nk\rj,\)i2e-6)] 



A comparison with ( p4| ) shows that 
Vs = V*- 



(48) 



(49) 



Note that we could have derived the slow-roll spectrum of 
from the exact spectrum ( ^ ) right from the beginning, 
by approximating it in the slow-roll regime. However, we 
have chosen to take the Bessel function/horizon crossing 
approach, because it is this approach which has been 
discussed in the literature. Let us note that the transfer 
function for C is unity. This means that the spectrum 
of ^ during the matter dominated era is identical to the 
spectrum at the end of inflation (region III) . 

We are mostly interested in the spectrum of the metric 
potential since this quantity appears in the calculations 
of the multipole moments, see Eq. (p^). If we assume that 
the Universe is matter dominated at the surface of last 
scattering, then the conservation law provides us with the 
relation ( = (5/3)$. Then, the spectrum of the Bardeen 
potential follows from (^ ) as: 



1-Ae + 26, 



(50) 



%^ [1 - 2e - 2(C + lnfcoh,|)(2e - 6)] . (51) 
257re 



These expressions are consistent with (4.3) and (5.1) 
of ijl^. The amplitude of scalar perturbations blows 
up when the slow-roll approximation becomes accurate, 
i.e. when e goes to zero. 

To end this section, let us make a last comment. It is 
clear from the previous considerations that we need the 
slow-roll approximation in region II only. In particular, 
this scheme of approximation is not needed in region III 



since the "exact" solution is known. However, one may 
wish to use it in region III also. Then, in this region, the 
Bardeen potential is given by $ ~ (Cs/2)e(l — 3e -|- 26). 
The long-wavelength transfer function, which allows to 
pass from the end of inflation to the matter dominated 
epoch can be expressed as r$ « [9/(25e^)](l -t- 6e — 46). 
Using the two previous formula, one can show that one 
recovers the spectrum given in Eqs. (|5^) and (|5|). How- 
ever, in principle, this method is not appropriate since 
we use an approximated solution whereas an exact one is 
available. 



B. Gravitational waves 

For gravitational waves, the same lines of reasoning 
can be applied. In region II, the effective potential can 
be written as: 



Ut{v) = a'H^ (2 - e) , 
and gives in the slow-roll limit 



3e 



(52) 



(53) 



Therefore the matching of sub- and superhorizon solution 
is again reduced to power-law inflation. The solution of 
/iT is similar to the one given in Eq. (|39|), where the 
effective index of the Bessel function is now given by: 



,(sr) _ 



(54) 



This solution can be used to find the constant Ct • Then, 
the power spectrum of gravitational waves reads 



(1 - 2Ce- 2elnfc) 



from which we deduce that: 



,(sr) 



-2e, 



,16 



aP ^ll,Hr-[l-2{C 



1 + lnA:o|?7*l)f 



(55) 

(56) 
(57) 



We see that there exists a crucial difference between den- 
sity perturbations and gravitational waves. In the case of 
gravitational waves, the ambiguity related to the choice 
of the matching time is not present. 

The amplitudes of scalar and tensor modes versus the 
slow-roll parameter e are displayed in Fig. ^ for 6 — e and 
in Fig. for 5 — 2e a.t leading and next-to-leading order. 
The first case is an approximation to the exact power-law 
result, the case 5 = 2e is the slow-roll approximation to 
Starobinsky's exact solution. 
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0.15 



FIG. 6. The scalar and tensor amplitudes from the 
slow-roll approximation for e — S. The scalar amplitude di- 
verges in the de Sitter limit e ^ 0. The leading order is drawn 
by full lines, the next-to-leading order by dashed lines. We 
have set A:o|?7*| = 1- 
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FIG. 7. The same as Fig. |, but for 5 = 2e. 



C. Multipole moments 

Let us first start with the calculation of Cf . We write 
the scalar multipoles as 



As 



(58) 



which defines the function gi{ns), cf. Eq. (p^). To com- 
pute cf at the next-to-leading order in the slow-roll pa- 
rameters we need to expand gi{ns) to first order around 
ns = 1, since ng itself is a linear function of e and S: 



ffi("s) .9i("S = 1) + ("S - 1)t^('^s = 1) 

dns 



2n 



■[l + (ns-l)(A-lnA:onss)] 



(59) 
(60) 



and \E'(a;) = dlnr(x)/da;. For the quadrupole we have 
L>2 « 1T463 and for large I, Di = 1 + ln(//2) + 0(l/l), 
due to ^'(/) = \nl + 0{l/l). Using Eqs. @) and (||), we 
find the scalar multipoles at next-to-leading order as 



0;2 TJ2 



1 



25e + 



l-2e-2(A +C)(2e-5) 



(62) 



In this equation the Doppler effect, the integrated Sachs- 
Wolfe effect, and the evolution of the transfer function 
are neglected. In the next section we will argue that this 
does not prevent the estimation of the slow-roll error. 

We now calculate R in the slow-roll regime at the lead- 
ing and next-to-leading order. The scalar quadrupole fol- 
lows from m2i) and reads: 



;2 tt2 

_ 'pi^* 
^ 75e 



l-2e-2(i:>2 + C)(2e-(5) 



+2(2e-5)lnp^ 

\v*\. 



where Da C « 0.4167. 

Let us now compute Cj. Using Eqs. (|lj) and 
find at next-to-leading order 



(63) 



we 



l-2(B + C + l-ln^ 



(64) 



where the number B is defined by 
B = dkk-^\n{k)I^{k)/ J^^ dkk-^I^ik) w 1.2878 so 
that B + C 0.5582. For higher tensor multipoles the 
numerical values of the constants in (|6^) are modified, 
but not the functional dependence on the slow-roll pa- 
rameters. In Fig. ^ the scalar and tensor quadrupoles at 
leading and next-to-leading orders are displayed for the 



case e 



6. 
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where 



Di = 1 -ln2 + «'(0 



ljrl/2 
l{l + l) 



FIG. 8. The scalar and tensor quadrupole moments from 
the slow-roll approximation for e = S. The leading order is 
(61) drawn by full lines, the next-to-leading order by dashed lines. 
We have set riss = |77»|. 
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Taking into account the expressions for C| and Cj 
given previously, we finally find the following expression 
for R in the slow-roll regime: 



R= 13.86e 



1 + 0.5504e - 0.8334(5 - 



(65) 



At leading order we recover the so-called consistency con- 
dition for slow-roll inflation which reads 



R = -6.93nT 



(66) 



This equation cannot be generalized by the use of ( pSj ) 
to a next-to-leading order equation, because it would in- 
volve the knowledge of the order O(e^) terms in nx- As 
discussed above, terms of that order are not meaningful 
in the slow-roll approximation. 

In Fig. H, the ratio R is displayed at leading and next- 
to-leading order for the two cases e — 6 and 2e = S. 




FIG. 9. The tensor to scalar ratio at leading order (full 
line) and at next-to-leading order for S = e (long dashed line) 
and 5 — 2e (short dashed line). The leading order is indepen- 
dent of S. We have set riss = |. 



V. DISCUSSION OF ERRORS 

The aim of this section is to quantify the magni- 
tude of the error introduced by the slow-roll approx- 
imation. For this purpose, we compare the slow-roll 
predictions with the exact results of power-law infla- 
tion. We explicitly test the following quantities: Q G 
{ns — 1, riT, As, At, Cf , CJ, i?}, i.e. quantities related to 
the power spectra and the quadrupole moments. 

We denote by Q the exact result of power-law infla- 
tion and by Q*-"-*, Q*-^-* the slow-roll results at leading 
and next-to-leading orders, respectively. The error is es- 
timated by calculating: 



=( ■0 - 



Q« - Q 



X 100% 



(67) 



Let us start with an estimate of the errors in the predic- 
tion of the spectral indices ns and tit. For the leading 



order slow-roll approximation m 



the error is e. 



(0) 



ns — 1 = and thus 

-ns-i ~ ''"T — 100%, except for de Sitter 
inflation. It is absolutely compulsory to use the next-to- 
leading order result for the spectral indices. We express 
the error as a function of 7. The best slow-roll approx- 
imation to a power-law model is given by e = S = j, 



and therefore n, 



(sr) 



log - 1 for this case. Thus from 
Eqs. (^, (|2^), (^, and (|5|) the error in the spectral 
indices from the slow-roll approximation is 



7 X 100% 



(68) 



Thus the next-to-leading order slow-roll approximation 
predicts the spectral indices with an error less than 1%, 
if 7 < 0.01 or 0.98 < ns < 1. 



Error in % 




FIG. 10. The error in the scalar quantities. The full 
lines are the quadrupole moments, the dashed lines are the 
amplitudes. The thin lines are the leading order corrections, 
the thick lines are the next-to-leading order corrections. 




0.05 0.1 0.15 

7 = e = (5 

FIG. 11. The same as Fig. but for the tensor amplitude 
and quadrupole moment. For the amplitude the error is the 
same as for the scalar sector, because e = 5. 

So far, except in some of the figures, we did not specify 
the pivot scale kg. We now choose to fix fco = a^^H^, i.e. 
fco is the mode that crosses the horizon at the time ry*, 
which is the time when we determine the values of the 
slow-roll parameters, where we fix them once and forever. 
It is easy to show that this amounts to take |r;^| ~ k^^ 
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in the argument of the logarithm which shows up in the 
equations of the previous section. StiU 77* remains to 
be fixed, which can be done most conveniently by fixing 
fcoTiss. In the following we will show that this choice is 
of physical relevance, because the accuracy of the slow- 
roll approximation can be improved by a clever choice 
of the pivot scale. We will discuss two cases. The usual 
convention is to choose fcoHss = 1- This corresponds to 
^phys _ ft,/ (6000 Mpc) today. Below we show that this 
leads to huge errors form the slow-roll approximation. 

To improve the precision of the slow-roll approximation 
we suggest to minimize the error in the region of the first 
acoustic peak, i.e. around I ~ 200. For this purpose we 
choose a pivot scale such that fcoHss ~ lOOe, which cor- 
responds to a physical wavenumber kf^' = /i/(22 Mpc) 
today. Let us start by analyzing the errors for the pivot 

The errors in the amplitudes and quadrupoles for the 
case of power-law inflation, e — S, are displayed in 
Figs. |l^ and ^ From these two plots, we can draw 
three conclusions. The first conclusion is that the error 
in the quadrupoles is larger than the error in the am- 
plitudes. This confirms the results already obtained in 
Ref. . The second conclusion is that it is not possible 
to obtain an error at the 1% level with the leading order, 
except for very small values of the slow-roll parameters. 
The third conclusion is that the accuracy of the next- 
to- leading order for the quadrupoles is better than 1% 
if 7 < 0.07, which corresponds to 0.85 < ns. Since the 
slow-roll approximation is more accurate for power-law 
model, it is reasonable to expect larger errors for more 
realistic models. 



100 




500 1000 



FIG. 12. The error from the next-to-leading order 
slow-roll approximation in the scalar multipoles Cf versus 
I for various values of the spectral index ns. The approxima- 
tion is best close to todays horizon since, in this figure, we 
have taken fcoriaa — 1, the most common choice. 

Let us now turn to the slow-roll errors in higher scalar 
multipole moments. Although we cannot obtain the ex- 
act result for the high-? multipole moments without mak- 
ing use of a Boltzmann code, we can nevertheless esti- 
mate the errors from the slow-roll approximation in this 
regime. As already discussed around Fig. |l| computing 



the scalar multipoles from Eq. ( [T^ ) is a bad approxima- 
tion for high I. Despite this fact, it is clear that for 
a given cosmological model the transfer function is the 
same for a power-law model and its slow-roll approxi- 
mation. Thus the only difference between the power- 
law and slow-roll multipole moments can arise from the 
convolution of this transfer function with different ini- 
tial spectra. We expect that this difference is small. To 
put it differently, C/(sr, t) — C;(sr, a) is large, whereas 
C/(sr, t) — C;{pl, t) ~ Ci(sr, a) — C;(pl, a), where t (a) 
denotes the use of the true (approximated) transfer func- 
tion and sr (pi) denotes the initial spectrum. Thus we 
use Eqs. (p^), (pO|), and (^) to obtain the errors for the 
scalar multipoles as a function of I, (|67|), which are dis- 
played in Fig. [ij. It shows that these errors are large 
and increase with I and \ns — 1|. 
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FIG. 13. The same as figure |l^, but with fcoHss = lOOe. 
This optimizes the accuracy in the region of the first acoustic 
peak. 

The reason for the large errors in the multipole mo- 
ments are the large errors in the spectral indices. This 
can be understood from the relations between the errors: 



100 
100 




(69) 



(70) 



where the prime denotes a derivative with respect to no,. 
The signs in front of the errors in the amplitude and in 
the spectral index are model dependent. For b = t the 
error in the amplitude has always positive sign, the error 
in the spectral index has a negative sign. For small values 
of I we may expand 

5;(ns) = gi{\) + g[{\){ns - 1) + 0[(ns - 1)'] (71) 

in Eqs. ( |69| ) and (|70|). Keeping only terms linear in ns — 1 
and terms linear in the errors we find 
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+ (A - lnA:onss)(l - "■s)e„. 



(i) 



-1 ' 



(72) 



where Di has been introduced in Eq. (|6l|). The error in 
the quadrupole moment is now easily understood from 
the last equation. As claimed above, the large error in 
the spectral index is responsible for the large error in the 
quadrupole moment. The contribution from the error in 
the spectral index always dominates. It is obvious from 
Fig. |l^ that the error increases with I and with |ng — 1|. 
Equation ( [72| ) underestimates the error for large I and 
|ns — 1|, due to the breakdown of the expansion ([7l|). 

These large errors, displayed in Fig. Q can be shifted 
to different multipoles by a change of the pivot scale fco- 
Therefore, one may hide part of the error from the slow- 
roll approximation in the cosmic variance. Inspection of 
Eq. ( [7^ ) suggests that the error is minimized for a given 
multipole if Di — ln(fco''iss)- For large / 



2 riss 



phys 



(73) 



where Hq is todays Hubble rate. We would like to know 
the position and hight of the first acoustic peak very accu- 
rate, which suggests to choose ^minimal error ~ 200. This 

gives fco = l/22/i-i Mpc « 1/31 Mpc for h = 0.71, the 
HST key project final value. Figure |l^ shows the er- 
rors at the next-to-leading order for various values of the 
scalar index with the new choice for the pivot scale. It 
can be seen clearly that the errors are highly suppressed 
around I « 200, as expected, but increase at lower and 
higher multipoles. The tiny bump between I « 200 and 
I « 500 is due to the fact that we plot the absolute value, 
in this region the error changes its sign. The new choice 
of the pivot scale allows to predict the multipoles in the 
range 2<l< 2000 for ng = 0.9 better than 10%, which 
was not possible with the pivot scale chosen previously. 
Nevertheless, the precision is not good enough to reach 
the 1% accuracy level (the error is 2.4% at I = 2000). 
In order to do so it is necessary to have ng > 0.93 or 
7 < 0.032. 




FIG. 14. The error in the tensor to scalar ratio in leading 
(thin line) and next-to-leading order (thick line). 

The error in the T/S ratio is displayed in Fig. |l^. In 
our special situation the pivot scale does not enter R, 



because e = S [see Eq. (|65|)]. We see that the error in 
R is less important than for the amplitudes and/or the 
quadrupoles. Therefore this suggests to use R to test 
the single scalar ficld/slow-roll paradigm. However, it is 
clear that any violation of the consistency check by the 
forthcoming data should be interpreted as a failure of 
this paradigm but not as the failure of inflation itself. In 
a more general situation where e S, the choice of fco 
does affect the error in R. Since we do not have an exact 
solution for scalar and tensors modes such that e S at 
our disposal, it is difficult to predict the corresponding 
effect. 

The errors from the slow-roll approximation displayed 
in Figs. |l^, |l^, 12, and |l^ should be compared with the 
cosmic variance. The cosmic variance is the variance of 
the best unbiased estimator for the multipole moments 

5Bcst(G) = 1/(2^+1) Er' 



25F 



where we 



have expanded the temperature fluctuations over the ba- 
sis of spherical harmonics, A(e) = X^zm o,imYim{e)- This 
expression is valid only if the afm's satisfy a Gaussian 
or a mildly non-Gaussian statistics. The corresponding 
error can be written as 



(cv) 



^^2pl(£il X 100% = J-^ X 100% . (74) 
Ci V 2/ + 1 



Over the whole range of the spectrum that will be mea- 
sured by a mission like Planck the cosmic variance is 
larger than 2%(^ k, 2000) and is « 7% at the first acous- 
tic peak (/ w 200). 

It is possible to reduce the cosmic variance by binning 
several multipoles together at the expense of decreasing 
the precision on the location of the multipoles. Therefore, 
we define an averaged multipole, C;, on the range [I — 
L,l + L] by 



C, 



1 



l+L 



2L + 1 



(75) 



j=l~L 



The central value I of each interval must be separated by 
2L -I- 1. For an incomplete sky coverage L is restricted 
by the form of the basis used to expand A(e). Below 
we do not take this issue into account and assume that 
the full sky is covered. In order to calculate the cosmic 
variance associated with the binned multipole, we define 
the estimator 



1 



l+L 



jU + 1) 



2L 



1 ^ 2j + l 



J 



(76) 



This estimator is clearly unbiased, {£[Ci)) = Ci, and it 
is very likely that this is also the best one although a 
rigorous proof is not presented here. Its variance can be 
expressed as 



1 



{2L 



l+L 

1)2 E 

' j=l-L 



2j + l 



(77) 
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Using that I {I 
arrive at 

(cv) ^ <^£{Ci) , 



l)Ci « const, if L is not too big, we as the cosmic variance for large 



since e 



2L + 1 



\ 



l+L 

E 

j=l-L 



100% w -J—. V — ^ X 100% 



2j + l 



(78) 

For L = this error reduces to the known expression for 
the cosmic variance (|7^) . In this case it is independent of 
the spectral index, whereas for L ^ this is true as fong 
as l{l + l)Ci is approximately constant within the range 
[I — L,l + L]. With the same approximation we find that 



the slow-roll error in C/ is just e^'' 



The cosmic variance for different multipoles and bin- 
ning is displayed in table For comparison we give the 
errors in the multipole moments from the slow-roll ap- 
proximation at leading and next-to-leading order for the 
pivot scale corresponding to the present horizon and we 
also present the errors at the next-to-leading order for the 
pivot scale corresponding to the scale of the first acoustic 
peak. We present the results for two values of the spec- 
tral index, corresponding to l(2)cr errors in the COBE 
measurement. 

For fco — ^i^s^i the error from the slow-roll approxima- 
tion in leading order dominates over the cosmic variance 
already at ^ = 10 for ng — 0.9 and L = 0. For the next-to- 
leading order and ns = 0.9 the slow-roll error dominates 
over the cosmic variance at / = 100, for any binning of 
multipoles. In the case of ns = 0.6 only the error in the 
quadrupole from the slow-roll approximation at next-to- 
leading order is smaller than the cosmic variance. Only 
for |ng — 1| < 10^^ the error from the next-to-leading 
order is below 1%, which is the order of magnitude of 
the cosmic variance in the Silk damping regime of the 
spectrum. 

For fco = l/22/i~^Mpc, the situation is improved. The 
slow-roll error and the cosmic variance are of the same 
order of magnitude up to Z = 1000 (ns = 0.9). However, 
this is not sufficient to hide the slow-roll error in the cos- 
mic variance for the whole spectrum. For ns — 0.6, the 
slow-roll error exceeds the cosmic variance at small and 
large scales and is hidden in the cosmic variance only for 
a narrow range of scales around fco ■ Let us point out that 
values of the tilt of the order of 0.1 and larger are a real- 
istic possiblity, as is clearly seen from the maximum like- 
lyhood fits to the recent BOOMERanG and MAXIMA- 1 
data. 

We conclude that for a general model of inflation only 
numerical mode-by-mode integration can presently pro- 
vide predictions for the CMBR with less than 1% er- 
ror, unless all slow-roll parameters are less than 10^^. 
To give an example let us consider the case of chaotic 
inflation with a potential V (p^, for which the slow- 
roll parameters are e « p/200,S k, [p — 2)/200, and 

^ « (p - 2)(p - 4)/(200p), giving ^ « 1 - (p + 2)/100. 
Thus for p — 2(4) the errors from the slow-roll approx- 
imation at next-to-leading order are of the same order 



0.01(0.02) 

and 5 « 0(0.01), which corresponds to rig"'' w 0.96(0.94). 
However, already for p = 6, the slow-roll approximation 
in the next-to-leading order leads to errors that exceed 
the cosmic variance at high multipoles (e ~ 0.03, 5 ~ 

0.02, C « 0.007, thus n^"'^ « 0.92). 

Forthcoming high precision missions, especially the 
MAP and Planck satellites, will only be limited by the 
cosmic variance up to Z « 1000 and I « 2000 respectively. 
Therefore predictions from inflationary models should be 
made such that the slow-roll error does not exceed the 
cosmic variance. We have shown that there are slow-roll 
models which cannot meet this requirement. 

Another implication is that the large errors in the pre- 
dicted multipoles render all attempts to reconstruct the 
inflationary potential difficult. The reason for this is that 
reconstruction usually assumes that the primordial spec- 
trum, instead of the multipoles, is measured to a high 
precision. We have shown in this work, that the errors 
in the prediction of the multipoles are easily an order of 
magnitude larger. A first attempt to go directly from 
inflation to the calculation of the multipole moments has 
been put forward by Grivell and Liddle recently. In 
our opinion a purely numerical approach to this funda- 
mental issue is not fully satisfactory — better analytic 
methods are needed. 





L 


I 






c 


impr 


0.9 





2 


63% 


15% 


0.4% 


8.1% 






10 


31% 


34% 


3.4% 


4.1% 






100 


10% 


68% 


15% 


0.35% 






1000 


3.2% 




39% 


1.0% 




2 


10 


14% 


34% 


3.4% 


4.1% 






100 


4.5% 


68% 


15% 


0.35% 






1000 


1.4% 




39% 


1.0% 




4 


10 


11% 


34% 


3.4% 


4.1% 






100 


3.3% 


68% 


15% 


0.35% 






1000 


1.1% 




39% 


1.0% 


0.6 





2 


63% 


73% 


8.6% 


56% 






10 


31% 




92% 


36% 






100 


10% 






4.6% 






1000 


3.2% 






22% 




2 


10 


14% 




92% 


36% 






100 


4.5% 






4.6% 






1000 


1.4% 






22% 




4 


10 


11% 




92% 


36% 






100 


3.3% 






4.6% 






1000 


1.1% 






22% 



TABLE I. Comparison of errors from the cosmic 
variance and the slow-roll approximation in leading and 
next-to-leading order with fcof"iss = 1 and in next-to- leading 
order with fcoJ^iss = lOOe, for various values of the spectral 
index ns, the bin width L, and the multipole I. The "— " 
denotes errors that exceed 100%. 
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